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For a general class of analytic f{R, RapR"^ , Rap-ysR"^''^) we discuss the gravitational lensing in 
the Newtonian Limit of theory. From the properties of Gauss Bonnet invariant it is successful to 
consider only one curvature invariants between the Ricci and Riemann tensor. Then we analyze 
the dynamics of photon embedded in a gravitational field of a generic f{R, RapR"''^ )-Gravity. The 
metric is time independent and spherically symmetric. The metric potentials are Schwarzschild-like, 
but there are two additional Yukawa terms linked to derivatives of / with respect to two curvature 
invariants. Considering first the case of a point-like lens, and after the one of a generic matter 
. . , distribution of lens, we study the deflection angle and the angular position of images. Though 

T-H . the additional Yukawa terms in the gravitational potential modifies dynamics with respect to the 

' General Relativity, the geodesic trajectory of photon is unaffected by the modification if we consider 

I only /(i?)-Gravity. While we find different results (defiection angle smaller than one of General 

, Relativity) only thank to introduction of a generic function of Ricci tensor square. Finally we 

^ ' can affirm the lensing phenomena for all /(i?)-Gravities are equal to the ones known of General 

' Relativity. We conclude the paper showing and comparing the deflection angle and position of 

, images for f{R, RapR"^)-Gr'Ayity with respect to the ones of General Relativity. 

^ ' PACS numbers: 04.25. Nx; 04.50. Kd; 04.80. Cc 
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I. INTRODUCTION 

^ ! 

Despite of all nice results of General Relativity (GR), the study ofpossible modifications of Einstein's theory of 
Gravity has a long history which reaches back to the early 1920s [lH6[. While the proposed early amendments of 
Einstein's theory were aimed toward the unification of Gravity with the other interactions of Physics, the recent 
^ ■ interest in such modifications comes from cosmological observations (for a comprehensive review, see Q)- In fact the 
T-H \ presence of the Big Bang singularity, the fiatness and horizon problems Q led to the statement that Cosmological 
. Standard Model, based on GR and Standard Model of Particle Physics, is inadequate to describe the Universe at 
' extreme regimes. Besides from Quantun Field Theory point view, GR is a classical theory which does not work as a 
I fundamental theory, when one wants to achieve a full quantum description of spacetime (and then of gravity). 
00 The astrophysical and cosmological observations usually lead to the introduction of additional ad-hoc concepts like 

Dark Energy/Matter if interpreted within Einstein's theory. In fact the principal physical aspects are the cosmic 
' acceleration and the fiat galactic rotation curves. These aspects could be interpreted as a first signal of a breakdown 
. . \ of GR at astrophysical and cosmological scales [9l-fl7| and led to the proposal of several alternative modifications of 
K*" ■ the underlying gravity theory (see [I^ for the review). 

While it is very natural (from the theoretical point of view) to extend Einstein's Gravity to theories with additional 
geometric degrees of freedom, (see for example [l9l - [2l| for general surveys on this subject as well as [2^] for a list of 
d ' works in a cosmological context), recent attempts focused on the old idea of modifying the gravitational Lagrangian 
in a purely metric framework, leading to higher order field equations. As such an approach is the so-called Extended 
Theories of Gravity (ETG) which have become a sort of paradigm in the study of gravitational interaction. They are 
based on corrections and enlargements of the Einstein theory. The paradigm consists, essentially, in adding higher 
order curvature invariants and minimally or non-minimally coupled scalar fields into dynamics which come out from 
the effective action of quantum gravity 23]. A sub class of ETG are the Fourth Order Gravities (FOG) where we do 
not consider further scalar fields but only the curvature invariants. 

The motivation to modify the GR come from the issue of a full recovering of the Mach principle which leads to assume 
a varying gravitational coupling. The principle states that the local inertial frame is determined by some average of the 
motion of distant astronomical objects [l^] • This fact implies that the gravitational coupling can be scale-dependent 
and related to some scalar field. As a consequence, the concept of "inertia" and the Equivalence Principle have to be 
revised. For example, the Brans-Dicke theory [25[ is a serious attempt to define an alternative theory to the Einstein 
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gravity: it takes into account a variable Newton gravitational coupling, whose dynamics is governed by a scalar field 
non- minimally coupled to the geometry. In such a way, Mach's principle is better implemented |25l-j27l|. As already 
mentioned, corrections to the gravitational Lagrangian were already studied by several authors [3,UlH|soon after the 
GR was proposed. Developments in the 1960s and 1970s (28l - [3^ . partially motivated by the quantization schemes 
proposed at that time, made clear that theories containing only a Ricci scalar square term in the Lagrangian were 
not viable with respect to their weak field behavior. Another concern which comes with generic higher order gravity 
(HOG) theories is linked to the initial value problem. It is unclear if the prolongation of standard methods can be 
used in order to tackle this problem in every theory. Hence it is doubtful that the Cauchy problem could be properly 
addressed in the near future, for example within the theories with inverse of Ricci scalar, if one takes into account 
the results already obtained in fourth order theories stemming from a quadratic Lagrangian (33l . |33 | . 

Besides, every unification scheme as Superstrings, Supergravity or Grand Unified Theories, takes into account 
effective actions where non-minimal couplings to the geometry or higher order terms in the curvature invariants are 
present. Such contributions are due to one-loop or higher loop corrections in the high curvature regimes near the 
full (not yet available) quantum gravity regime [ISl- Specifically, this scheme was adopted in order to deal with 
the quantization on curved spacetimes and the result was that the interactions among quantum scalar fields and 
background geometry or the gravitational self-interactions yield corrective terms in the Hilbert-Einstein Lagrangian 

From a conceptual viewpoint, there are no a priori reasons to restrict the gravitational Lagrangian to the linear 
function of the Ricci scalar, minimally coupled with matter [37| . Since all curvature invariants are at least second 
order differential, the corrective terms in the field equations will be always at least fourth order. That is why generally 
we call them higher order terms (with respect to the terms of GR). 

The idea to extend Einstein's theory of gravitation is fruitful and economic also with respect to several attempts 
which try to solve problems by adding new and, most of times, unjustified ingredients in order to give a self-consistent 
picture of dynamics. The today observed accelerated expansion of the Hubble flow and the missing matter of as- 
trophysical large scale structures, are primarily enclosed in these considerations. Both the issues could be solved 
changing the gravitational sector, i.e. the l.h.s. of field equations. The philosophy is alternative to add new cosmic 
fluids (new components in the r.h.s. of field equations) which should give rise to clustered structures (dark matter) 
or to accelerated dynamics (dark energy) thanks to exotic equations of state. In particular, relaxing the hypothesis 
that gravitational Lagrangian has to be a linear function of the Ricci curvature scalar, like in the Hilbert-Einstein 
formulation, one can take into account an effective action where the gravitational Lagrangian includes other scalar 
invariants. In summary, the general features of ETGs are that the Einstein field equations result to be modified in two 
senses: i) geometry can be non-minimally coupled to some scalar field, and / or ii) higher than second order derivative 
terms in the metric come out. In the former case, we generically deal with scalar-tensor theories of gravity; in the 
latter, we deal with higher order theories. However combinations of non-minimally coupled and higher-order terms 
can emerge as contributions in effective Lagrangians. In this case, we deal with higher-order-scalar-tensor theories of 
Gravity. 

Due to the increased complexity of the field equations in this framework, the main amount of works dealt with some 
formally equivalent theories, in which a reduction of the order of the field equations was achieved by considering the 
metric and the connection as independent fields (36l - l40| . In addition, many authors exploited the formal relationship 
to scalar-tensor theories to make some statements about the weak field regime, which was already worked out for 
scalar-tensor theories more than ten years ago [4l|. Moreover other authors exploited a systematic analysis of such 
theories were performed at short scale and in the low energy limit [i^ - lisl . [sol Isil [ssl |59| . In particular the FOG has 
been studied in the so-called Newtonian Limit (case of weak field and small velocity) and in the Weak Field Limit 
(case of Gravitational waves) (60j . In the first case we found a modification of gravitational potential, while in the 
second one the massive propagation of waves. After these preliminary studies it needs to check the new theories by 
applying them in the realistic model. Then the galactic rotation curve [13, [6ll - l63j or the motion of body in the Solar 
System have been evaluated. 

It is worth remembering that one of the first experimental confirmations of Einsteinian theory of Gravity was the 
deflection of light. Since then, gravitational lensing (GL) has become one of the most useful successes of GR and 
it represents nowadays a powerful tool [64]. In this paper, we investigate the equation for the photon deflection 
considering the Newtonian Limit of a general class of f{X, Y, Z)-Gravity where / is a unspecific function oi X = R 
(Ricci scalar), Y = RajjR'^^ (Ricci tensor square) and Z = R^p-ysR'^^^^ (Rieman tensor square). 

The starting point of paper is the metric tensor solution of modified Einstein equation [5l| and we analyze the 
dynamics of photon. Particularly in Section II we introduce the /(X, Z)-Gravity and its field equations resuming 
shortly the Newtonian limit approach, while the Section III-A is devoted to study of photon traveling embedded in a 
gravitational field. In this section we approach first the deflection angle by a point-like source and after we reformulate 
the lensing problem by a generic matter distribution (Section III-B). Moreover we rewrite the equation lens (Section 
III-C) and find the corrections to the position of images in the case of point-like lens. Finally in Section IV there are 
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comments and conclusions. 



II. THE NEWTONIAN LIMIT OF FOURTH ORDER GRAVITY 

Let us start with a general class of FOG given by the action 



A = / (Tx^/^ 



fiX,Y,Z) + XCr^ 



(1) 



where / is an unspecified function of curvature invariants. The term Cm is the minimally coupled ordinary matter 
contribution. In the metric approach, the field equations are obtained by varying ((!]) with respect to g^j^. We get 



/ 



2fzR^o.pjR.''^^ -^fzR^^^J-cp = XT^, (2) 



where T,, 



1 Sj^C^) 



is the energy-momentum tensor of matter, fx 



dX 



fy 



df. 

dY 



fz 



dZ ' ^ 



and 



X — 8ttG^. The convention for Ricci's tensor is = R'^ ^av, while for the Riemann tensor is R°' p^i, = r|^^ + .... 
The affinities are the usual Christoffel symbols of the metric: F^^ — \g^'^ {gaa.p + ff^o-.a — dap. a)- The adopted 

signature is (H ) (for details, see [11]). 

The paradigm of the Newtonian limit starts from the development of the metric tensor (and of all additional 
quantities in the theory) with respect to the dimensionless quantity v but considering only first term of tt- and 
ij-component of metric tensor g^i, (for details, see [lOl). The develop of the metric is the following 



ds^ = (1 + 2^)df - (1 - 2'^)5^jdx'dx' 



(3) 



where $ and ^ are proportional to . The set of coordinates^ adopted is — (/:, x^, x^, x'^). The curvature 
invariants X , Z become 



z - + . . . 



(4) 



and the function / and its partial derivatives (fx, fxx, fY and fz) can be substituted by their corresponding Taylor 
develop. In the case of / we have 



fix, Y, Z) ^ fiO) + /x(0)X(2) ^ l/xx(0)X(2)' + /;f (0)XW + /y(o)r(4) + /z(0)zw + 

and analogous relations for derivatives are obtained. 
From the lowest order of field equations ([2|) we have 



(5) 



/(O) - 



(6) 



while in the Newtonian Limit (oc w^) we have^ 



^ Here we use the convention c = 1. 

^ The Greek index runs from to 3; the Latin index runs from 1 to 3. 

^ Throughout the paper we assume always fx (0) > 0, and therefore we may set fx (0) = 1 without loss of generality. 
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(A-TO2^)A$ 



7712 



mi +2m2 



3m 



(7) 



A* = /rfVg2(x,x')(^^- "^6m'"^ A..jX(^)(xO 

^ (A-mi2)x(2) = mi^A-p 

where X^^'' is the Ricci scalar at Newtonian order, p is the matter density and Q2 is the Green function of field 
operator A — m2^. The quantities m^^ are linked to derivatives of / with respect to the curvature invariants X, Y 
and Z 



mi 



7712 



'3/xx(0)+2/y(0)+2/z(0) 



(8) 



/i'(0)+4/z(0) 

By solving the field equations ([T]), if > for i =^ 1,2, the proper time interval, generated by a point-like 
source with mass M, is (for details, see (sol. Isij) 



1 1 e^^il^l 4 e^^^l''! 



1 le^^il'^l 2e"^2l''l 



5,.jdx'dx^ (9) 



where = 2GM is the Schwarzschild radius and /i; = -^/jroi^. The field equations ([7]) are valid for any values 
of quantities to^^, while the Green functions of field operator A — rrii'^ admit two different behaviors if m^^ > or 
rrii'^ < 0. The possible choices of Green function, for spherically symmetric systems {i.e. tj,i(x,x') = Gi{\x — x'|)), 
are the following 



e.(x,x') = 



1 e-^-l"-" ' 
in |x— x'l 

1 cos /jj |x— x' I +sin /j^ |x— x' I 
47r |x— x'l 



if TTli^ > 

if m,^ < 



(10) 



The first choice in pop corresponds to Yukawa-like behavior, while the second one to the oscillating case. Both 
expressions are a generalization of the usual gravitational potential (cx Ixj"^), and when rrii'^ — 00 (i.e. 
fxx{0) ,/i'(0) ,/z(0) — > from the ([8])) we recover the field equations of GR. Independently of algebraic sign of 
rrii^ we can introduce two scale lengths fii~^. We note that in the case of / (X)-Gravity we obtained only one scale 
length (/Ui~^ with /y(0) = /z(0) = 0) on the which the Ricci scalar evolves [HO, [Ell, but in /(X, F, Z)-Gravity we 
have an additional scale length /i2~^ on the which the Ricci tensor evolves. 

Often for spherically symmetric problems it is convenient rewriting the metric ([9]) in the so-called standard coordi- 
nates system^ (the usual form in which we write the Schwarzschild solution) . By introducing a new radial coordinate 
f = Ixl as follows 



the relativistic invariant © becomes^ 



1 - 2*(r) 



(11) 



-Ml'' 



dt-" 



pir + 1 



2(M2r+l) 



dr^-r^dn (12) 



where d^ = d9^ + sin 6 dcf)^ is the solid angle and renamed the radial coordinate f. 



* Generally the set of coordinates (t, r, d, 4>) are called standard coordinates if the metric is expressed as ds^ = gtt(t, f) dt^ +grr(t, r)dr^ — 
r^dQ while if one has ds^ = gtt{t, x) dt^ + gij {t, x)dx^dx^ (like the solution ^) the set {t, , , x^) is called isotropic coordinates [57| . 
^ The metrics II12I I and ^ represent the same space-time at first order of rg/r. 
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III. THE GRAVITATIONAL LENSING BY f{X, Y, Z)-GRAVITY 
A. Point-Like Source 

The Lagrangian of photon in the gravitational field with metric (fT2|) is 



C 



1 _ !± E{r) ]p- [1 + ^ A(r) | r" - r 
r J \ r 



sin^ 06^ 



(13) 



where S(r) = 1 + S(r) = 1 + 1 e^^i'' - | e-^2^ A(r) =^ 1 + A(r) = 1 - g-'^i'' - ^iifsphli g-z^^r and the dot 

represents the derivatives with respect to the afHne parameter A. Since the variable 9 does not have dynamics (9 — 0) 
we can choose for simplicity 9 — 'k/2. By applying the Euler-Lagrangian equation to Lagrangian (|13p for the cyclic 
variables t, 6 we find two motion constants 



dt 



and respect to A we find the "energy" of Lagrangian^ 



(14) 



= -J 



£ = 

By inserting the equations (fll)) into we find a differential equation for f 



(15) 





1 J2" 


l + ^A(r) 





(16) 



r+ is the solution for leaving photon, while f_ is one for incoming photon. Let tq be a minimal distance from the 
lens center (Fig. [3]). We must impose the condition r±(ro) = from the which we find 



J2 = 



1-^S(.„) 

Now the deflection angle a (Fig. [3]) is defined by following relation 



(17) 



-dA = 



■ dr 



■ dr = -TT + 2 



ro ' + 



■ dr 



Ao 



)(iA 



)dA 



(18) 



where Ao is the value of A corresponding to the minimal value (ro) of radial coordinate r. By putting the expressions 
of J, and r+ into ((T8)) we get the deflection angle 



dr 



■ A(r 



I{r) ro-' 



(19) 



The I I13I I is a quadratic form, so it corresponds to its Hamiltonian. 



6 



which in the case Vg/r ^ l'' becomes 



(20) 



where 



-^A'i,M2(''o) = / , dr (21) 

From the definition of S and A we note that in the case /(AT, F, Z) — > AT we obtain J^^i^. ^i^ifo) — > 0. In a such way 
we extended and contemporarily recovered the outcome of GR. 

The analytical dependence of function ^;ii,^2(''o) from the parameters /ii and 112 is given by evaluating the integral 
(I2ip . A such as integral is not easily valuable from the analytical point of view. However this aspect is not fundamental, 
since we can numerically appreciate the deviation from the outcome of GR. In fact in Fig. [T] we show the plot of 
deflection angle ((20|) by f{X, Y, Z)-Gravity for a given set of values for fii and /i2- The spatial behavior of a is ever 
the same if we do not modify 112- This outcome is really a surprise: by the numerical evaluation of the function 
-^/Ji- ^2 (''o) one notes that the dependence of fii is only formal. If we solve analytically the integral we must find a /ii 
independent function. However, this statement should not be justified only by numerical evaluation but it needs an 
analytical proof. For these reasons in the next section we reformulate the theory of GL generated by a generic matter 
distribution and demonstrate that for /(X)-Gravity one has the same outcome of GR. 




0.0^ — ^ — ^ — ^ — < — ^ — ^ — ^ — < — ^ — ^ — ^ — < — ^ — ^ — ^ — < — ^ — ^ — ^ — u 
2 4 6 8 10 

'0 



FIG. 1: Comparison between the deflection angle of GR (solid line) and one of f{X, Y, Z)-Gravity (dashed line) (|20|) for a fixed 
value /X2 = 2 and any /xi. 



B. Extended Matter Source 



In this section we want to recast the framework of GL for a generic matter source distribution p(x) so the photon 
can undergo many deviations. In this case we leave the hypothesis that the flight of photon belongs always to the 
same plane, but we consider only the deflection angle as the angle between the directions of incoming and leaving 
photon. Finally we find the generalization of GL in /(X, y, Z)-Gravity including the previous outcome of deflecting 
point-like source (and resolving the integral (|21I) V 

The relativistic invariant ([3]) is yet valid since we consider the superposition of point-like solutions. Indeed we can 
generalize the metric (O by the following substitution 



^ We do not consider the GL generated by a black hole. 
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1 + i e^'^il'^^'^'l — I e^'^^lx-x'i 



X — i g-Pilx-x'l _ 2 g-/j2|x-x'| 



(22) 



This approach is correct only in the Newtonian hmit since a such Umit correspond also to the linearized version of 
theory. Obviously the f{X, Y, Z)-Gravity (like GR) is not linear, then we should have to solve the field equations 
with a given p. 

By introducing the four velocity = = {vP ,vl) the flight of photon, from the metric Q, is regulated by the 
condition 



gapu w 



(l + 2$)u"^-(l-2«')|u|2 = 



then u'^ is given by 



(23) 



1 - 2* 



u ,u 



1 + 2$ 

In the Newtonian limit we find that the geodesic motion equation becomes 



(24) 



|upV($ + f) - 2uV* • u = 



and by supposing |up — 1 we can recast the equation in a more known aspect 



(25) 



u = -2 



Vi* + -V($-*) 



(26) 



where V_l 



1^ • V 1 1^ is the two dimensions nabla operator orthogonal to direction of vector u. In GR we 



would had only u = — 2 Vj^$ since we have VE* = $. In fact the field equations ([7]) are corrects [51] if we satisfy a 
constraint condition among the metric potentials $, ^' as follows 



A($- *) = 



nil 



1712 



dVg2(x,x')Ax'X(2)(x') 



(27) 



We can affirm, then, that only in GR the metric potentials are equals (or more generally their difference must be 
proportional to function |x|~^). The constraint ((27|) has been found also many times in the context of cosmological 
perturbation theory 52- 56]. 

The deflection angle ([T8| is now defined by equation 



a 



(28) 



where and A / are the initial and final value of affine parameter |66| . For a generic matter distribution we can not a 
priori claim that the deflection angle belongs to lens plane (as point-like source), but we can only link the deflection 
angle to the difference between the initial and final velocity u. So we only analyze the directions of photon before and 
after the interaction with the gravitational mass. Then the ([25]) is placed by assuming a = Au = — u^. From 
the geodesic equation ([^5]) the deflection angle becomes 



V_l1' + -V($-1') 



(29) 
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The formula ([29)) represents the generaUzation of deflection angle in the framework of GR. By considering the photon 
incoming along the z-axes we can set = (0,0,1). Moreover we decompose the general vector xeR'^ in two 
components: ^eK^andzeR. The differential operator now can be decomposed as follows V = Vj^ + z9z = V^+zc^z, 



while the modulus of distance is |x — x'| = y |C ^ ^'P + ^ ^0^ ^(C; C ^ ^')- Since the potentials $ , ^' <^ 1, 
around the lens, the solution of ([26l) with the initial condition = (0, 0, 1) can be expressed as follows 



u = (Cl($,^'),e)($,*),l + e)($,*)) 
and we can substitute the integration with respect to the afflne parameter A with z. In fact wc note 



(30) 



dX 



dz 



dz 



dz/dX l + 0($,*) 



dz 



and the deflection angle becomes 



a 



V,-($ + 4') + zaz($-*) 



dz 



From the expression of potentials ([2^ we find the relations 



(31) 



(32) 



$ + * = -2 G f d2|*'dz'-4t^l- + 2 G f d^Cdz \^ii^'"> , e"^^^(«^«'-^'-') 

^ -z.z') ^ A( f .f' .z.z') 



$ - $ = - ^ / d^^'dz' P'^"'-"') 



g-Ml A(|',|'',z,z') _ g-/i2A(|',|'',z,z') 



and the equation ((^ becomes 



(33) 



a =2Gl\^edz'dz(^^^^-^'^ 



A(e,e,z,z') ^ A(C,e,z,z') 



3 A. A(f,e,z,z') 



1 + Ail A(C, I*', z, z') e-^^^(«'«''^'^') - 1 + M2A(C, C', z, z') e-^^^(€.€^ 



(34) 



,2,2: ) 



In the case of hypothesis of thin lens belonging to plane (x, y) we can consider a weak dependence of modulus 
A(^, z, z') into variable z' so there is only a trivial error if we set z' = 0. With this hypothesis the integral into z' 
is incorporated by definition of two dimensional mass density S(^') — J dz' p{^\ z'). Since we are interesting only to 
the GL performed by one lens we can extend the integration range of z between (— cxd, oo). Now the deflection angle 
is the following 



4G / d2^'S(r 



(35) 



where 



' ^ Jl + /.2A(e,$^z,0)) ^_^^^(,-,,-.,,,o^ 
A(C,e,^,0)' 



(36) 



The last integral in p4p is vanishing because the integrating function is odd with respect to variable z. The expression 
(|35p is the generalization of outcome (j20l) and mainly we found a correction term depending only on the fi2 parameter. 
In the case of point-like source S(C') = M5'^^^(^') we find 
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a 



2r„ 



1 

I\ 



1^1 -^^.(e.o) 



1 



(37) 



and in the case of f{X,Y,Z) f{X) (i.e. ^2 oo and ^i_12{£,t^') ~> 0) we recover the outcome of GR a — 
2 r-g ^/If p. From the theory of GL in GR we know that the deflection angle 2 rg/ro is formally equal to 2 rg/\^\ if we 
suppose To = |^|. Besides both ro, |^| are not practically measurable, while it is possible to measure the so-called 
impact parameter b (see Fig. [3]). But only in the first approximation these three quantities are equal. 

In fact when the photon is far from the gravitational source we can parameterize the trajectory as follows 




r = Vt^ + b^ 



arctan : 



(38) 



and from the definition of angular momentum (1141) in the case oi t ^ b we have 



b/t^ 



1 + b^/f^ 

By using the condition (jl7p r±(ro) = we find the relation among b and rg^ 



(39) 



b = 



ro T 



l-^^2(ro) 



(40) 



justifying then the position ro = |^| in the limit rg/r <C 1 (but also Vg/ro <C 1). 

In Fig. [2] we report the plot of deflection angle ([37|. The behaviors shown in figure are parameterized only by ^2 
and we note an equal behavior shown in Fig. [TJ With the expression ([37]) we have the analytical proof of statement at 
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FIG. 2: Comparison between the deflection angle of GR (solid line) and of /(X, Y, Z)-Gravity (dashed line) (|37p for 0.2 < 
M2 < 2. 



the end of previous section. In fact in the equation (\37^ we have not any information about the correction induced in 
the action ([1} by a generic function of Ricci scalar {fxx 7^ 0). This result is very important if we consider only the 
class of theories /(X)-Gravity. In this case, since fi2 ^ 00, we found the same outcome of GR. From the behavior in 



The constant T is dimensionless if we consider that A is the length of trajectory of photon. In this case without losing the generality 
we can choose T = 1. 
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Fig. [2] we note that the correction to outcome of GR is deeply different for tq — )■ 0, while for rp — )■ oo the behavior 
([57]) approaches one of GR, but the deviations are smaller. This difference is given by the repulsive correction to the 
gravitational potential (see metric ([9])) induced by f{Y, Z). Only by leaving the thin lens hypothesis (the lens does not 
belong to plane z = 0) we can have the deflection angle depending by /ii ([M)) . In fact in this case the third integral 
in (p4)) is not zero. Then in the case of thin lens we have a complete degeneracy of outcomes in the /(X)-Gravity: 
all f{X)- Gravities are equivalent to the GR. If we want to find some differences we must to include the contributions 
generated by the Ricci tensor square. But also in this case we do not have the right behavior: the deflection angle is 
smaller than one of GR: f{X,Y, Z)-Gravity does not mimic the Dark Matter component if we assume the thin lens 
hypothesis. 



C. Lens equation 

To demonstrate the effect of a deflecting mass we show in Fig. [3] the simplest GL configuration. A point-like mass is 
located at distance Dql from the observer O. The source is at distance Dqs from the observer, and its true angular 
separation from the lens L is /3, the separation which would be observed in the absence of lensing (rg = 0). The 
photon which passes the lens at distance ro ~ 6 is defiected with an angle a. 

Since the deflection angle ([50]) is equal to (1571) , for sake of simplicity we will use the " vectorial" expression. Then 
the expression ([571) . by considering the relation (|in)) . becomes 



a = 2rg 



(41) 



The condition that this photon reach the observer is obtained from the geometry of Fig. [3l In fact we find 



f3 = e-^a (42) 
Uos 

Here D^s is the distance of the source from the lens. In the simple case with a Euclidean background metric here, 
Dls = Dos — DoL] however, since the GL occurs in the Universe on large scale, one must use a cosmological model 
[66j . Denoting the angular separation between the deflecting mass and the deflected photon as 9 = b/DoL the lens 
equation for f{X, Y, Z, )-Gravity is the following 



^J-{9)]e^ - 13 9 -0E^ = (43) 

•2rgDLS 



where 9e — \/ j^J^ j,^^ is the Einstein angle and 



Ti9)= I d,(i±ifH££4^^e--^^--^) (44) 

Jo ^{9^ + 7jy 

Since we have < 9'^F{9) < 1 (Fig. |4]) we can find a perturbative solution of (j43|) by starting from one in GR, 
gGR ^ -l}±y/fs^+ieE l^ j^^^ assuming 9 ^ 9%^ + 9* and neglecting 9*^F{9*) in (gS]) we find 

e = 9rT , , H0r)0r" (45) 

and in the case of /3 = we find the modification to the Einstein ring 



9 = ±9f 



■J'{9e) 



(46) 



In Fig. [5] we show the angular position of images with respect to the Einstein ring. Both the deflection angle and 
the position of images assume a smaller value than ones of GR. Then the corrections to the GR quantities are found 
only for the introduction in the action ([T]) of curvature invariants Y (or Z), while there are no modifications induced 
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FIG. 3: The gravitational lensing geometric for a point-like source lens L at distance Dql from observer O. A source S at 
distance Dos from O has angular position /3 from the lens. A light ray (dashed line) from 5* which passes the lens at minimal 
distance ro is deflected by a\ the observer sees an image of the source at angular position 9 = b/DoL where b is the impact 
factor. Dls is the distance lens - source. 
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FIG. 4: Plot of function e^T{e) for 1 < fi2DoL < 10. 

by adding a generic function of Ricci scalar X. The algebraic signs of terms concerning the parameter /i2 are ever 
different with respect to the terms of GR in ([9]) and they can be interpreted as a " repulsive force" giving us a minor 
curvature of photon. The correction terms concerning the parameter /ii have opposite algebraic sign in the metric 
component gtt and gij ^ and we lose their information in the deflection angle (j32p . 

In both approaches we find the same outcomes /.ti-independent because the matter source (in our case it is a point- 
like mass) is symmetric with respect to z-axes and we neglect the second integral in Obviously for a generic 
matter distribution the deflection angle is defined by and the choice of second derivative of function of Ricci 
scalar is not arbitrary anymore. 
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FIG. 5: Plot of the Einstein ring (solid line) and its modification (|46p in the f{X, Y, Z)-Gravity for 1 < fi2DoL < 10 (dashed 
line). 

IV. CONCLUSIONS 

In this paper we compute the study of GL when a FOG is considered. Among the several theories of fourth order 
we consider a generic function of Ricci scalar, Ricci and Riemann tensor, but by using Gauss-Bonnet invariant it is 
adequate to consider only a f{X, F)-Gravity. 

We start from the outcome of previous paper [5l| about the point-like solutions in the so-called Newtonian limit 
of theory and formulate the deflection angle and the angular positions of image. The study has been evaluated on 
two steps: in the first one we consider a point-like source and by analyzing the properties of Lagrangian of photon we 
obtain a correction to the outcome of GR depending apparently on two free parameters of theory. But by plotting, 
only numerically, the new angular behavior ()20p with respect to the minimal distance ro we note that the correction 
term does not depend on the parameter fii (Fig. [T]). In the second step we start by more general geodesic motion 
and we reformulate the deflection angle for a generic matter distribution. In the case of an axially symmetric matter 
density we obtain the usual relation between the deflection angle and the orthogonal gradient of metric potentials 
([35| . Otherwise we flnd that the angle is depending also onto the travel direction of photon (|34|) . Particularly if there 
is a z-symmetry the deflection angle does not depend explicitly on the parameter /ii but we have only the correction 
term induced by ^2 ([36)1 . 

From the deflnition of /ii and 112 ^ we note that the presence of function of Ricci scalar {fxx{0) 7^ 0) is only 
in 111. Then if we consider only the f{X)-Gravity (fj,2 -^00) the geodesic trajectory of photon is unaffected by the 
modification in the Hilbert- Einstein action. Our analysis is compatible with respect to the principal outcome shown 
in the paper [g^I . Instead if we want to have the corrections to GR it needs to add a generic function of Ricci tensor 
square into Hilbert- Einstein action. But in this case we flnd the deflection angle smaller than one of GR (f20|) or p7|) . 
Obviously the same situation is present also in the Einstein ring (|46p , where the new angle is ever lower than the one 
of GR (Fig. [5]). The mathematical motivation is a consequence of algebraic signs of terms containing the parameter 
/i2 in the metric ([9]) . In fact they are ever different with respect to the terms of GR in ([9]) and they can be interpreted 
as a "repulsive force" giving us a minor curvature of photon trajectory, instead the correction terms containing the 
parameter /xi have opposite algebraic sign in the metric components gtt and gij ([9]) and we lose their information in 
the deflection angle (|32|) . 

A similar outcome has been found for the galactic rotation curve (s^l, where the contribution of f{Y) in the action 
gives us a lower rotation velocity proflle than the one of GR, but with a no trivial difference. In fact in galactic 
dynamics we are studying the motion of massive particles and in this case we flnd the corrections induced also by 
/(X)-Gravity. Then if we can estimate the weight of the corrections (induced by /(X)-Gravity) to the Ricci scalar 
for the galactic motion, from the point of view of GL we have a perfect agreement with the GR. Only by adding 
f{Y) in the action we induce the modiflcations in both two frameworks, but we do not flnd the hoped behavior: the 
flat galactic rotation curve and a more strong deflection angle of photon. Also for a photon bending we need a Dark 
Matter component. Moreover if we consider f(X,Y)- Gravity for the GL we need a bigger amount of Dark Matter 
then in GR. A such conclusion qualitatively does not differ by the one about the galactic rotation curves [63| . 

Also after the analysis of GL we can affirm, as for the galactic rotation curves, it remains a hard challenge to 
interpret the Dark Matter effects as a single geometric phenomenon. 
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